boundary_types = p,p
boundary_values = 0,7
initial_gradients =      -0.2, -0.5, -0.7,
reflection_path = (1,2),(2,2),(0,3),(2,2),(1,3)
grid_size = 5
max_radius = 70%
